#model
model{
   #likelihood
   for(i in 1:nGirls){ Girls[i] ~ dnorm(muGirls, precGirls)}
   for(i in 1:nBoys){ Boys[i] ~ dnorm(muBoys, precBoys)}


   sGirls <- sqrt(1/precGirls)
   sBoys <- sqrt(1/precBoys)
   diff <- muBoys-muGirls
   
   #priors
   muGirls~dnorm(50, .00000000000000000001)
   muBoys~dnorm(50, .00000000000000000001)
   precGirls~dgamma(.001,.001)
   precBoys~dgamma(.001,.001)
}



#data
list(nGirls=26, nBoys=24,Girls=c(
 72, 69, 80, 90, 77, 84, 59, 72, 62, 82, 66, 72, 96, 80, 74, 81, 82, 81, 69,
 68, 85, 69,106, 95, 60, 64),
Boys=c(
 86, 79, 83, 80, 78, 64, 89, 89, 81, 86, 75, 71,122,104,103, 65, 92,108,
 86, 85, 88, 78, 91,103))

#inits
list(muGirls=50,muBoys=50,precGirls=1,precBoys=1)


#Second approach

#model
model{
      #likelihood
      for(i in 1:n){Peabody[i]~dnorm(mu[i],precision)
      mu[i] <- a + b*Sex[i] }

      #priors
      a~dnorm(0,.00000000000000001)
      b~dnorm(0,.00000000000000001)
      precision~dgamma(.001,.001)

      sd <- 1/sqrt(precision)
}


#data
list(n=50, Peabody=c(
 72, 69, 80, 90, 77, 84, 59, 72, 62, 82, 66, 72, 96, 80, 74, 81, 82, 81, 69,
 68, 85, 69,106, 95, 60, 64, 86, 79, 83, 80, 78, 64, 89, 89, 81, 86, 75, 71,
 122,104,103, 65, 92,108, 86, 85, 88, 78, 91,103),Sex=c(0,0,0,0,0,0,0,0,0,0,0,
 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))

#inits
list(a=0,b=1,precision=.5)